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Abstract 



In this letter we develop a method to construct all the integrals of motion 

>' 
lO ' of the SU(p) Haldane-Shastry model of spins, equally spaced around a circle, 

interacting through a 1/r 2 exchange interaction. These integrals of motion 

■*-)- ' respect the Yangian symmetry algebra of the Hamiltonian. 
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The past few years have seen an extensive study of exactly solvable quantum many body 
systems with l/r 2 -interactions. The simplest member of this family is the SU(p) Haldane- 
Shastry model (HSM) with Hamiltonian |[]||: 

describing N particles with an internal spin degree of freedom that can take on p different 
values, residing on equally spaced sites on a ring: {zj = exp (^p j)}- For p — 2 it describes 
spin | particles. The operator Py permutes the spin of two particles at sites i and j. The 
energy levels of this model turn out to have huge degeneracies (beyond the regular global 
SU(2) symmetry), signaling the presence of a large non-trivial symmetry algebra. In Ref. 
H this algebra was identified as the Yangian, a Hopf algebra introduced by Drinfel'd in 
1986 ||. It describes the elementary excitations of this model: spinous. These spinons obey 
semion fractional statistics. 

The fact that this Yangian algebra commutes with the Hamiltonian hints at the in- 
tegrability of this model. However, the traditional method of proving integrability, i.e. 
construction of a set of commuting extensive hermitean operators {Hi,H 2 , . . .}, so-called 
invariants, has so far been unsuccessful. Invariants up to H4 have been 'guessed' |||3|]. Mi- 
nahan and Fowler introduced a set of invariants that commute with the Hamiltonian, 
based on operators introduced by Polychronakos 0. However, the generating function for 
this set is essentially the trace of the transfer matrix and thus contains only elements of the 
Yangian algebra.Q 

In this letter we will construct a set {H n } of extensive operators that commute among 
themselves and with the Yangian. In order to do this we have to consider a more general 
dynamical model in which the particles are allowed to move along the ring: the Calogero 



x The authors of Ref. || claim to have found the Hamiltonian H2 in their third order invariant, 
but in reproducing their algebra we did not find any such term; in fact we only recovered yangian 
operators. 



Sutherland model (CSM) with an internal degree of freedom. This model, which has been 
studied in J7||| has the following Hamiltonian: 

When A — ► oo or equivalently ft — > 0, the particles 'freeze' into their classical equilibrium 
positions, and barring some subtleties we recover the spin Hamiltonian Hi. 

The reason for this diversion through the dynamical model to obtain the constants of 
motion is the following: the so-called quantum determinant of the transfer matrix, an object 
that commutes with the Yangian algebra and therefore a natural candidate for the generating 
function of the constants of motion, happens to be scalar in the spin model (i.e. when h — ► 0), 
as we shall see below. But in the general dynamical model, this is not the case, and by 
carefully taking the limit A — ► 00 we can isolate a generating function for the {H n }. 

Let us first review the role of the Yangian algebra in the dynamical model. A more 
extensive treatment can be found in P,[T0|j. The integrability of the CSM is based on the 
existence of the transfer-matrix T al '(u) that commutes with the Hamiltonian: 



00 \ 
T ab (u) = 5 ab + Y —rTZ b 

n=0 " 

N 

Lij = 8ijZjd z . + (1 - Sij)X6ijPij, (3) 

where Xj b , a,b = l,...,p acts as \a)(b\ on the spin of particle i, and 9ij = z z J; z . . This 
transfer matrix satisfies the Yang-Baxter equation: 

Ro ,(u - v)T°(u)T°'(v) = T°'(v)T°(u)R 00 ,(u - v), (4) 

with Rq i(u) = u + APoc an d T°(u) = T(u) ® 1 and T° (u) = 1 ® T(u). For the purposes 
of this letter we will deal with another form of the same transfer matrix. Introduce the 



following representation of the so-called DnnH-operators [jlO 



A = hzidg. + 7j = hzid Zt + -^ (wij - sgn(i - j))K i:j 

Z{ ~r Zj . . 

w i:j = _ (5) 

Z{ Zj 



Kij is the operator that permutes the spatial co-ordinates of particles i and j. These Dunkl- 
operators commute: 

[D i ,D j }=0, (6) 

but are not covariant under permutations: 

[K iti+1 , D k ] = tik^i,i + l 

Ki,i+iD{ — D i+ xKi t i + i = A, (7) 

defining a so-called degenerate affine Heche Algebra. In terms of these Dunkl operators we 
can define a transfer matrix that also obeys the Yang-Baxter equation: 

TV)=fl + ^^V--f 1 + ^^ L V (8) 

V u-Dj V u-D N ) 

It satisfies eq. (f|) trivially since the {Di} commute amongst themselves and commute with 

the Pqj, since the latter only act on spin degrees of freedom; furthermore ( 1 + _ % J is the 

elementary transfer matrix with spectral parameter u — Di. To retrieve T°(u) from T°{u) 

we have to apply a projection II to T° that replaces every occurrence of K i>i+ i by Pi <t+ \ 

once ordered to the right of an expression (this is equivalent to having the unprojected 

operator act on wavefunctions that are symmetric under simultaneous permutations of spin- 

and spatial co-ordinates) []XQf| . We will drop the "0" subscript on T(u) from here on. 

Normally the conserved quantities are derived from a Taylor expansion of the trace of the 

transfer matrix. In this case that just gives us combinations of Yangian operators, elements 

of the symmetry algebra. This set doesn't even contain the Hamiltonian. As pointed by 

various authors [II,3.IIJ, there is another quantity that commutes with the Hamiltonian, 



derivable from the transfer matrix: the quantum determinant, 

Det q (T(u))= J2 e{a)T lai {u-\(jp-l))T 2ff2 {u-\(p-2))-.-T p(Tp {u). (9) 

UdSp 



It satisfies [T(u),Det q (T(u))] = 0. It has been computed in [fjlj as: 

Det„(T(u)) = IT Det 9 (f («))n = II (^+M) n. (10) 



where (making the dependence on h explicit): 

A(u,h) = l[(u-D l (h)). (11) 

i 

Now obviously: 

[A(u,h),f(v,h)] = 0. (12) 

This holds since the D^s commute with each other and the Po/s . The projector doesn't 
get in the way since a product of projections is the projection of the product — both T(y) 
and A(u) are symmetric under simultaneous permutation of spin- and spatial co-ordinates 
PH| | . The eigenvalues of A(u) are also known: for every partition \n\ there is an eigenvalue: 

jM(u) = l[(u-hn j -\V --£-)) (13) 

3=1 / 

We notice that as h — > 0, i.e. in the limit of the HSM, all eigenvalues become identical 
and A(u, 0) is a multiple of the identity operator. Thus no non-trivial constants of motion 
are contained in A (it, 0). Nevertheless let us study fll2"D for small h. Writing T(v,h) = 
J2 n h n f n (v); A{u,fi) = ?2 n h n A n (u): 

0= [f(v,h),A(u,h)] 
= [f (v), A (u)] + h {[f (v), Ax(«)] + [fi(v), A («)]) + 0{h 2 ). (14) 

The 0(h ) term is trivially 0. The rest of this letter will focus on the vanishing of the 0(K) 
term. As we shall show below: [Tx(v ), A (m)] = 0. Therefore we have the important result: 
[T Q (v), Ai(it)] = 0, i.e. the 0(h) term in A(u, n) commutes with the transfer matrix (and 
therefore the Yangian) of the Haldane-Shastry spin model as well. Furthermore it will also 
become apparent that 

[A 1 (w),A 1 (u')] = 0. (15) 

So we can take Ai (it) to be the generating function of the constants of motion of the HSM! In 
order to establish these results we first need to prove the following corollary: Zid Zi A (u) = 
when evaluated with the particles at their equilibrium positions, i.e. Zj = exp (^pj)- 



From eq. (jTT|) we have A (m) = Y\i(u — jA. Since we know that A (m) is scalar we 
can evaluate it by having it act on any convenient state, e.g. the one where all particles 
have identical values for their internal degree of freedom (for p = 2 we would say: all spins 
pointing up). I.e. the permutations reduce to 1. This has been done in [ |10|j : 



A (m) = det(w - 6) 
\z~ 
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;i-w 



(16) 



Then using d x detA(x) = Tr[A 1 d x A(x)]det[A(x)] we have: 

1 



-a,A n (n)=Tr 



w-e 



d Zi Q 



det[u-6] 



E u~ {n+1) Ti [e n d Zi Q] det[u - 0]. 

n=0 



(17) 



Now evaluate the trace in a basis where is diagonal. This can clearly be done for h — > 
and Zj = exp (^pj)- has eigenvectors ip n , where (ip n )j = -4^ exp (^jn\ with eigenvalue 
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n). Then 
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Tr[0 m «9 2i 0] = Y.md z m n )\ m (^^ - n 



(18) 



Using that 
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we find: 
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Therefore <9^A(w) = at Zj = exp (^pj). 

From expanding A(u,h) to C?(7i) in ( |TT]) we have 



= 0. 



(20) 



N (i-1 \ IN 

aiw^e n(«-7i) ^ n ( m -7,) 

*=i Vi=i / v=*+ 1 



(21) 



Then, with the corollary and the fact that [7$, jj] = for all i,j (the Dunkl algebra (0) is 
satisfied for Ti = as well): 

N fi-1 \ / N 

[a 1 («),a («)]=e nc«-7i) [^ Zi ,A («)] n c«-7i: 

AT A-l \ /at 

= £ II(«-7i) {^Ao(«)} II(«-7*))=0 (22) 

With this result we find: 

1 1 . 1 * , . 1 _ 2 , 



h-. — Ai(u)^ — + 0(/r 



A(«,n) A («) A (w) A (n 

= XTT " K—^—tMu) + 0(h 2 ). (23) 

A (w) (A («)) 

This can be checked by multiplying the LHS and RHS by A(u, ft). Now let us expand T(v, ft) 
to 0(h): 

N ( \P \ 1 N 1 N 

t(v,h) = H l + ^r = ■ xr -rrII(«-A + APoi) = TT TlI(«-7i + APoi) + 

i= i \ v-DiJ A(v, ft) i=1 A (v) i= i 

+/ M-(xtt) Ai( v ) n (« - 7* + AP W ) 

[ \A (u)/ i= i 

' E ( S(f - 7i + APo,) I *0« [ fi (« ~ 7i + APo,) ] 1 + 0(ft 2 ) 



= f («) + KT r (v) + 0{% 2 ) (24) 

It is now obvious how [Ti(v) , A (u)] = 0. The contribution from the first term in curly 



brackets in (^4j) vanishes by virtue of eq. (22). As far as the second term is concerned: 



the 7j commute amongst each other and with the Pqj, and [zid Zi , A (w)] = as we showed 
before. 

So far we have established that Ai(u) respects the Yangian symmetry, but we also need 
to show that it is a good generator of constants of motion in that it commutes with itself 
at a different value of the parameter u: [Ai(w), Ai(tt')] = 0. It will be enough to prove this 
on the space spanned by the Yangian Highest Weight states (YHWS). These states are, as 
their name implies, the highest weight states of a representation of the Yangian algebra. 



All other states in the model are generated by acting on these YHWS with the elements of 
the Yangian algebra, i.e. the transfer matrix. Since Ai(u) commutes with T (v) for any u, 
[Ai(w), Ai(u')] = will therefore also hold on the non-highest weight states. First of all we 
should note that Ai(u) and Ai(u') leave the space of YHWS invariant. This follows from 
the fact that all such states |T) are annihilated by Tq(v) with a > b (see ref. [[TJJ). But 
since T^ 6 (i>) commutes with Ai(ii), T^ ( Ai(w)|r)J will also be for a > b. 

The proof that Ai(w)and Ai(tt') commute hinges on the existence of an operator that 
commutes with both these Ai's and is non- degenerate. Such an operator is Tq P (v). Its 
eigenvalues are given by ||10|| : 



P 1 (v + l)---P p _ 1 (v + l) 



with v = v/X (25) 



Pt(v) ■ ■ ■ Pp^iv) 

The polynomials Pi(v), i = 1 . . .p—1 characterize the representation of the Yangian. As was 
found in fL0"l , every degenerate supermultiplet in the SU(p) HSM (i.e. every representation 



of the Yangian) can be represented by a string of iV + 1 binary digits or 1, called a motif. 
The first and last (entry and N) are always 0. For an SU(p) model the string cannot 
contain more than p—1 consecutive ones. A set of k — 1 consecutive l's is called a k- 
string. To generate the polynomials P, replace very in the motif by ')('. We then have 
groups of l's enclosed in parentheses. The polynomial Pkiy) has zeroes at | + rightmost 
member of a fc-string. Let us consider an example to clarify this: for SU(3), N = 14 the 
motif 011010001100100 has four 1-strings, two 2-strings and two 3-strings. This implies 
polynomials P x {v) = (v - 5f)(t> - 6|)(U - 10|)(v - 13±), P 3 (v) = (v - 4|)(U - 12|) and 
P,(v) = {v-2\){v-9\). 

The eigenvalues ( p5| ) are obviously independent rational functions of v, and Tq P (v) is non- 
degenerate. If |T) is a YHWS with motif F then Ai(w)|r) and Ai(w')|r) both are scalar mul- 
tiples of |T) since they have the same Tq P (w) -eigenvalue ([Tq P (v), Ai(m)] = [Tq P (v ), Ai(w')] = 
0). So in this YHWS-space both Ai(u) and Ai(w') are diagonal and thus commute. 

In the remaining part we will reproduce the integrals of motion that have already been 
found f5|||, and point out some subtleties in their construction. As is customary for the 
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Heisenberg model with nearest neighbor exchange we take Ti(u) = 4- ln(Ai(w)) rather than 
Ai(tt) to be the generating function for the integrals of motion, so that the invariants will 
have an additive spectrum. When expanded in powers of u it reads: 

N l 

v^u) = nE — 7-^ -n to o(h) 

~{u- hzid Zi - ji 

oo ( N ra-1 

= E u- (n+1) nE E^m^r^n 

n=0 I i=l p=0 



= Y. u ~ {n+1)H m (26) 

n=0 

where we have reinserted the projection operator that turns Kij — > Pij, when ordered to the 
right of an expression. We have worked out the first few H's. With z^ = z% — Zj we have: 

N 

#i = E^ = p 



i,3 % 



ijfc z ij z jk z ki 1 ij 

ijk Z ij Z jk z ki 4 

^ = E'^^^(^-1)-2E' {^-)\p i3 -1) + ^±H 2 (27) 

ijfcj z ij z jk z kl z h ij \ z ij z jiJ 6 

The prime on the summation symbol indicates that the sum should be restricted to distinct 
summation-indices. To compute the previous expressions we normal ordered the Zid Zi to the 
right in eq. ( pop and then put Zj = exp \jffj)- The identity WijWjk + WjkWki + w^iWij = —1 
which lies at the heart of the integrability of these -^--models is very useful in the reduction 
of these expressions. P indicates the momentum and we will show its interpretation later on. 
Notice the absence of Yangian operators as well as terms containing both permutations and 



derivatives. The expressions in eq. fl27[ ) can be seen to coincide with those reported previously 
@,f] lending credibility to this way of deriving the integrals of motion. Unfortunately for 
large n it becomes prohibitively complicated to compute H n . 



2 We should point out a correction in eq. (7) of Ref. |J where — 4-^2 should be replaced by +|i?2- 



We have an alternative way to verify the validity of these constants of motion as well. 
We will proceed to compute the eigenvalues of the operators H n and compare these to the 
'rapidity' description of the eigenvalues in |J. We will constrain ourselves to the SU(2) case 
to simplify the algebra. 

As is well known |T2j the roots of A(u, Ti) — i.e. the poles of the transfer matrix T(u, K) , see 
eq. ([24]) — are given by the solutions of so-called Bethe Ansatz equations, which only depend 
on the two-particle phase-shift. In the case of the CSM it is 7rAsgn(&4 — k-z). Notice that 
this phase-shift only depends on the ordering of the momenta k\ and &2- This is why these 
models are interpreted to describe an ideal gas of particles with statistics that interpolates 
between bosons (A = 0) and fermions (A = 1). In the dynamical model (0) the particles have 
charge and spin. Therefore we get two coupled sets of 'Nested' Bethe Ansatz equations — for 
the general case p 7^ 2 there are p equations. They have been presented in |~3| : 



kiL = A <{ J2 nsgniki - kj) + - 



M 

lixli - ix Y^ sgn(/cj - A c 

a=l 



= (k$ + -tai) L (28) 

_ N 

n J2 s g n ( A « - A p) + 27rJ a = 7r^sgn(A a - fc). (29) 

We have reinstated L, the circumference of the circle to get the dimensions correctly. Notice 
how h drops out of these equations due to the fact that the full Hamiltonian is scale invariant 
(a peculiarity of the -^-type potentials). So rather than sending S -» we should let 
A — > 00. There are N equations defining the {hj} (one for every particle) with charge 
quantum numbers {/«}. Furthermore we have M equations defining the auxiliary momenta 
{A a }. M is the number of particles with a spin j. The {J a } are their spin quantum numbers. 
The I's and J's are distinct integers or half-odd integers depending on the parity of N and 
M-N. 



This changes the invariant in a harmless manner, but this is relevant for comparing the eigenvalues 
of the operators in that article and the ones that we will find later on. 
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Furthermore we should restrict the Hilbert space to states that only carry spin excita- 
tions, and no charge excitations (these elastic modes unfortunately don't acquire a gap as 
ft — » 0)[|. As in for instance the lD-Hubbard model we accomplish this by leaving the charge 
quantum numbers in their groundstate configuration and only exciting the {J a }- Let us 
therefore first analyze the absolute groundstate which has M = y (for N even), so there 
are twice as many fc's as A's. The fc's and J's are consecutive and spaced by one unit. 
Then eq. (|29"D tells us that between every two A's there must be two fc's. In a spin excited 
state we have M < y and by leaving openings in the J-distribution, we can have more 
than two fc's sit between every pair of A's. Notice that this equation doesn't fix the value 
of the A's, just their positions with respect to the fcj's. From eq. ( p8|) we learn that when 
we order the {fcj} such that fcj < kj for i < j: ki ~ A f i — — ^) = kf. There is however 
an O(t) = 0{h) correction through the A's. Whenever a A sits between two fc's they will 
be drawn together by t. This information is contained in 5k t . Now, in the same way that 
the constants of motion are contained in the (D(ft)-term in A(w,ft), their eigenvalues are 
determined by the O(j-) corrections to the fc's. As A(w, ft) has eigenvalues Y\f = i(u — ki), 
T(u, ft) = £ In A(u) =T + hT l + ... must have eigenvalues £i ^- = Ei u _ fe? ii,; fc . - So for 
ri(rx) acting on a state characterized by a set {A a } we find its eigenvalue: 

oo i N 

n=0 " i=l 

oo 1 

- E ^iMM- (so) 

n=0 " 

We will now label the A a 's by m a , their positions relative to the fc's, i.e. if A a sits between k r 
and fc r +i, then m a = r. We see that the m a have to be at least two units apart since there are 
at least two fc's between A's according to eq. (p9|). Now writing |sgn(A a — fcj) = ^(Aq, — fc») — ~ 
(6* is the step function), we have: 



3 This is why we don't use eq. (13) which gives a much more direct expression for the eigenvalues 
of A(u, ft); however, one doesn't know a priori whether an eigenvalue belongs to a pure spin- or 
charge excitation, or a mixture of both. 
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N N / 7V + In"' 1 f M M 

i=l t=l V Z 7 I Z a=l 



M 



const + 2J 



o=l 



. i=l 



= const + £ £m( m a) (31) 

o=l 

For small n we can evaluate e n (m a ) exactly: 

e x {m a ) = m a 

£2(^0) = m a (m a - N) 

, A^ 2 -l 
£3(^0) = 2 m a( m a - N)(2m a - N) ^ — m a 

( N 2 - 1 \ 
e 4 (m a ) = e 2 (m a ) I e 2 (m a ) H — 1 (32) 

These results coincide nicely with the numerical values of ||, when we interpret the m a as 
the rapidities of the HSM! We notice that it is indeed consistent to interpret the momentum 
term P in (|27|) in Hi and H 3 as J2a m ct, i>e. the degree of the polynomial wavefunction. 

In conclusion we have outlined a method to obtain the constants of motion of the HSM as 
a strong coupling limit of the CSM with particles with internal degrees of freedom. Although 
the task to actually obtain the invariants is quite cumbersome, it can be done in principle. 
Given the relatively simple structure of the invariants we expect there to be some technique 
that could simplify the computation considerably. The construction of integrals of motion 
presented in this letter, provides us with extensive operators that commute with each other 
and the Yangian symmetry algebra. By computing eigenvalues of the invariants through the 
Nested Bethe Ansatz, and comparing them with previous numerical results we provided 
evidence for the validity of the approach. It would be interesting to analyze the cause of the 
miraculous absence of terms containing mixtures of permutations and derivatives. 

This work was supported in part by NSF Grant No. DMR922407. 
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